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ABSTRACT 



a,: 

6 

Atomic hydrogen loss at the top of HD 209458b's atmosphere has been re- 
cently detected (Vidal-Madjar et al. 2003). We have developed a 1-dimensional 
model to study the chemistry in the upper atmosphere of this extrasolar "hot 
jupiter". The 3 most abundant elements (other than He), as well as 4 parent 
molecules are included in this model, viz., H, C, O, H 2 , CO, H 2 0, and CH4. The 
higher temperatures (~ 1000 K) and higher stellar irradiance (~6x 10 5 W m~ 2 ) 
strongly enhance and modify the chemical reaction rates in this atmosphere. Our 
two main results are that (a) the production of atomic hydrogen in the atmo- 
sphere is mainly driven by H2O photolysis and reaction of OH with H2, and is 
not sensitive to the exact abundances of CO, H2O, and CH4, and (b) H2O and 
CH 4 can be produced via the photolysis of CO followed by the reactions with H 2 . 



Subject headings: planetary systems — radiative transfer — stars: atmosphere — 
stars: individual (HD 209458) 
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1. INTRODUCTION 

Since the discovery of the first extrasolar planet, 51 Peg b, in 1995 (Mayor & Queloz 
1995), a total of 102 planets have so far been discovered (e.g., Butler et al. 2003; Udry et al. 
2002, and references therein) and analyzed statistically in order to characterize the formation 
environment (Fischer et al. 2002; Santos et al. 2003). The formation of gas giants is thought 
to be complete in 10 Myr, before the disappearance of the gaseous stellar accretion disk, at 
distances > 5 AU from the parent star. They are then pulled to their present positions by 
tidal interaction between the gas disk and planet (e.g., Pollack et al. 1996; Ward 1997). 

An edge-on planet provides a unique opportunity to investigate the planetary atmo- 
sphere. HD 209458b is such a planet, providing the first extrasolar planetary detection using 
the light curve obtained during a planetary transit of its parent star (Charbonneau et al. 
2000; Henry et al. 2000). The orbital parameters were accurately determined by Charbon- 
neau et al. (2000), Henry et al. (2000), and Mazeh et al. (2000). Strong absorption lines 
are required to make an atmospheric detection and Seager & Sasselov (2000) theoretically 
characterized the most prominent absorption features, viz., Na I and K I doublet resonance 
and He I 2 3 S-2 3 P triplets. Charbonneau et al. (2002) detected the Na I doublet at 589.3 nm 
in HD 209458b with ~ 4a confidence level. Following this, Vidal-Madjar et al. (2003) made 
the first observation of the extended upper atmosphere of HD 209458b with a ~ 4a detec- 
tion of the H I atomic hydrogen absorption of the stellar Lyman-a line. They reported an 
absorption of ~ 15 ± 4% and claimed this should be taking place beyond the Roche limit, 
thus implying hydrodynamic escape of hydrogen atoms from HD 209458b's atmosphere. 

The temperature and UV flux of close-in planets are high. This motivates us to study 
the chemistry that may be important in this "hot jupiter". In this paper, we consider a 
simple hydrocarbon/oxygen chemistry model to determine a source of atomic hydrogen in 
the atmosphere of HD 209458b and represents the first effort to investigate the UV enhanced 
chemical processes in "hot jupiters". 

2. MODEL 

Our model is based on the four parent molecules H 2 , CO, H 2 0, and CH 4 and is a 
derivative of the Caltech/JPL KINETICS model for the Jovian atmosphere. HD 209458b is 
orbiting at a distance of 0.05 AU. As HD 209458 is a GO solar-type dwarf star, it is justified 
to use the solar spectrum. We expect the atmosphere to have a temperature, ~ 1000 K, and 
UV flux (< 1800 A), ~ 2 x 10 15 photons cm 2 s 1 . By comparison the UV flux at Jupiter 
is a factor of ~ 10 4 lower. 
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Hydrocarbons The hydrocarbon photochemical scheme used here is a simplified version of 
the Jovian atmospheric model described in Strobel (1973) and Gladstone et al. (1996). The 
photodissociation of CH 4 and the subsequent reactions of the species with hydrogen produce 
all the other hydrocarbons present in a Jovian-type atmosphere. For lower temperatures and 
weaker stellar irradiation, the main source of H is from H 2 and CH4 photodissociation and 
the main sink via C2H2, which acts as a catalyst in recombining H. HD 209458b receives 
much greater stellar irradiation and is therefore much hotter than Jupiter. In this case, the 
formation of H is greatly enhanced by photolysis of H 2 and reactions between O and OH 
radicals and H 2 . The sink for H is more complex (see § 3.3). 

Oxygen O is similar in abundance to C and represents a cosmic abundance of these 
species. The atmospheric H 2 abundance will be controlled mainly by the comparative 
richness of these two species. The amount of H 2 and CO governs the amount of atomic 
oxygen present and the related reactions are of great interest. The oxygen related reactions 
are taken from Moses et al. (2000). 

Model atmosphere Three models are investigated in this paper. Our standard reference 
model using solar abundances is shown in Fig. 1 (Model A). We also consider 2 other cases, 
Model B and Model C, in which H 2 and CO abundances are 10 times lower, respectively. 
We have taken the 1 bar level to be "0" km, and all heights are referenced from this level. 
The temperature-pressure profile and chemical abundances are based on Seager et al. (2000). 
The temperature decreases from the bottom to the top of the atmosphere. The abundances 
of CO and H 2 are, 3.6 x 10~ 4 and 4.5 x 10~ 4 , respectively. These values are similar to 
solar abundances. The CH4 abundance is determined by the thermodynamic equilibrium 
chemistry in the deep atmosphere. We adopt the value 3.9 x 10~ 8 , which is the lower end of 
the model by Seager & Sasselov (2000). The temperature-pressure profiles are not certain, 
because global circulation and high temperature condensation (Seager & Sasselov 2000; 
Sudarsky et al. 2000) are not included in generating the model atmosphere. Nevertheless, 
the present standard reference model is accurate enough for a first order understanding of 
the chemistry and characterizes a source of H in the upper atmosphere. The eddy diffusion 
is proportional to n~ Q (n is number density) where a is taken to be ~ 0.6 — 0.7. 



3. RESULTS 

The principle results for our three models will be described in three sections. Section 
3.1 deals with OH and O radicals, followed by C0 2 and CH4 in section 3.2. The important 
question of hydrogen production in relation to H 2 is addressed in section 3.3. 
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Fig. 1. — Vertical profile of temperature, total density, and constituent number density. 
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3.1. OH and O radicals 

Fig. 2 shows the OH and O radicals in our models. OH and O are the most important 
radicals as they drive most of the chemical reactions. O is the most important element in 
facilitating the formation of OH radicals in the water-poor atmosphere. These species are 
produced photochemically, for example, 

CO + h// -> C + O, 

H 2 + his -»• H + OH. 

CO photolysis is an important source of O and H2O is the main source of OH. Though 
0( 1 D) is not as abundant as O (< 10~ 5 [O]), it can react with H 2 to produce a similar 
amount of OH. O is increasing with altitude as a consequence of H 2 and CO photolysis. 
OH is increasing with altitude until it starts decreasing at ~ 10 nbar. The decline of OH 
above 10 nbar is due to OH photodissociation. We see the mixing ratio of OH radicals is 
not sensitive to the abundance of CO and H 2 0. With an order of magnitude change in H 2 
(Model B) or CO (Model C), OH is changed only by a factor of < 3. However, O is sensitive 
to both CO and H 2 concentrations and preferentially forms OH. From Fig. 2 we see that 
OH is not sensitive to H 2 abundance. The abundance of H 2 depends on the comparative 
richness of cosmic C and O. Under high stellar UV irradiation, a fraction of CO will be 
photodissociated. The resulting O will react with H 2 to form OH, which eventually forms 
H 2 by reacting with H 2 . Therefore, the abundance of OH radicals is determined by the 
amount of O in the system. It does not matter whether O is in the CO or H 2 reservoir. 



3.2. C0 2 and CH 4 

C0 2 is formed via the reaction of CO and OH, 

OH + CO -> C0 2 + H 

Fig. 3 shows the vertical profiles of C0 2 for our three models. The C0 2 mixing ratio is 
enhanced in the upper atmosphere. At pressures of ~ 10 nbar for Model A, the C0 2 mixing 
ratio is only about 2 orders of magnitude less than its progenitor, CO. The C0 2 abundance 
in the model is rather insensitive to the abundance of H 2 0. An order of magnitude decrease 
in H 2 results in only a factor of ~ 3 decrease in C0 2 abundance (cf. Models A and B in 
Fig. 3). However, C0 2 abundance varies approximately linearly with the abundance of CO 
(cf. Models A and C in Fig. 3). 
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Mixing Ratio 



Fig. 2. — Comparison of mixing ratios of OH and O radicals in Models A, B, and C. 



-7- 



The formation of CH 4 is initiated by the downward flux of C atoms produced in the 
photolysis of CO in the upper atmosphere. This obtains the following sequence of reactions: 

CO + hz/^C + O, 

C + H 2 + M -> 3 CH 2 + M, 

2 3 CH 2 -> C 2 H 2 + 2H, 

C 2 H 2 + H + M -> C 2 H 3 + M, 

C 2 H3 + H 2 — > C 2 H 4 + H, 

C 2 H 4 + H + M -> C 2 H 5 + M, 

C 2 H 5 + H^2CH 3 , 

CH 3 + H + M — > CH 4 + M. 

Fig. 3 shows the vertical profiles of CH 4 in our models. We see that the CH 4 mixing ratio is 
increasing by a factor of 5-100 from the bottom to the 0.1 mbar level (cf. Fig. 3). Above this 
level, CH 4 rapidly decreases due to photodissociation. Some of the C is eventually converted 
to C0 2 , whose mixing ratio increases while that of CH 4 decreases. 

The CH 4 mixing ratio is increased by a factor ~ 2 when we lower the H 2 abundance 
by an order of magnitude (cf. Models A and B in Fig. 3). We suggest this increase is due 
to less UV shielding by water above. The CH 4 mixing ratio is decreased by an order of 
magnitude when we lower the CO abundance by an order of magnitude. The reason is that 
CO photolysis is the source of the C in CH 4 . 



3.3. H and H 2 

Fig. 3 shows the mixing ratios for H and H 2 0. The most striking features are the 
production of H (all three models) and the production of H 2 (Model B). In our 1-D model, 
the production rate of atomic hydrogen is not sensitive to the exact abundances of CO and 
H 2 0. With an order of magnitude change in the abundance of either CO or H 2 0, the atomic 
hydrogen changes by only a small factor, ~ 1 — 2. This implies the production of H in three 
models is limited by the availability of UV photons. The H production is also not sensitive 
to the abundance of CH 4 . CH 4 abundance has been increased to be as high as CO, and 
the H mixing ratio is only changed by a small factor. A more comprehensive discussion will 
be given in a separate paper (Liang et al. 2003, in preparation). The mixing ratio of H 
exceeds 10% at the top of atmosphere. At the top of the atmosphere (< 1 nbar), H 2 will 
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be photolyzed and is a source of H. This atomic hydrogen will fuel the hydrodynamic loss 
process, as observed by Vidal-Madjar et al. (2003) and is discussed in more detailed in the 
companion paper by Parkinson et al. (2003, in preparation). 

In a water-poor atmosphere (e.g., Model B), CO will be driving the photochemical 
reactions to form H 2 0: 

CO + h/y -> C + O, 
O + H 2 -»• OH + H, 
OH + H 2 ^H 2 + H. 

H 2 and H are the net products. The produced H 2 will be recycled and is an important 
source of OH radicals and H atoms. Near the top of atmosphere, a large fraction of H 2 is 
destroyed due to the high UV bombardment. However, CO is more stable. 

Fig. 4 shows the production rate of H and the photolysis rate of H 2 0. The rates are not 
sensitive to the abundance of CO, but are sensitive to the abundance of H 2 0. By comparing 
the H 2 photolysis rate with H production rate, it is evident that production of H is mainly 
driven by H 2 photolysis and the reaction of OH with H 2 . Below ~ 1 mbar, HCO plays a 
role in the removal of H via 

CO + H + M^HCO + M, 
HCO + H^CO + H 2 . 

In the upper atmosphere, H atom recombination and reactions with CH and 3 CH 2 will drive 
the loss of H via 

2H + M -> H 2 + M, 

and 

H + CH^C + H 2 , 
H + 3 CH 2 -> CH + H 2 , 
C + H 2 + M -> 3 CH 2 + M. 



4. CONCLUSION 

We have considered a series of possible chemical reactions using various models for a 
"hot jupiter". We have shown the mechanism for producing the atomic hydrogen. The 
production of H is not sensitive to the abundances of CO, H 2 0, and CH 4 . Lowering H 2 
or CO an order of magnitude changes the concentration of H by only a factor of < 2. 
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Mixing ratio 



Fig. 3. — Comparison of mixing ratios of H, CH 4 , H 2 0, CO, and C0 2 in Models A, B, and 
C. 
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Fig. 4. — Production rate of H (dark lines) and photolysis rate of H2O (gray lines) in Models 
A, B, and C. 
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However, the production rate of H is sensitive to the temperature profile. A 30% change in 
the temperature will result in ~ 50% change in the H concentration. 

Our calculations show that the H mixing ratio at ~ 1 mbar is ~ f -3 and exceeds 10% in 
the top of the atmosphere. Being less gravitationally bound, the atomic hydrogen formed at 
the top of atmosphere can escape hydrodynamically as putatively suggested by observations 
of Vidal-Madjar et al. (2003). Since these close-in gas-rich giant planets are probably tidally 
locked, circulation may be important in transporting heat because the temperature gradient 
can be as high as 1000 K over the entire planet, which in turn implies wind speeds of a few 
km s _1 (Showman & Guillot 2002; Cho et al. 2003). Therefore, it is of interest to simulate 
the differences in chemical processes between the day and night sides and the global transport 
of heat and mass as well. 

We thank Dr. R. L. Shia and Dr. M. Gerstell for helpful comments and discussions. The 
support of NASA Grant NAG5-6263 to the California Institute of Technology is gratefully 
acknowledged. This work represents one phase of research carried out at the Jet Propulsion 
Laboratory, California Institute of Technology, under contract to the National Aeronautics 
and Space Administrations. 
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